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ABSTRACT 


A general  formalism  is  developed  to  describe  the  nonlinear  evolution 
associated  with  the  parametric  decay  of  an  intense,  coherent  electromagnetic 
wave  into  an  electrostatic  wave,  its  second  harmonic,  and  scattered  electro- 
magnetic waves  in  a homogeneous  plasma.  The  effects  of  pi^>  depletion  are 
neglected  and  it  is  assumed  that  all  waves  are  coherent.  Two  classes  of 
solutions  aze  found.  One  class  is  explosively  unstable,  while  the  other  con- 
sists of  growing  aperiodic  oscillations.  The  evolution  of  modulational  modes 
is  discussed. 
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I.  INTRODUCTION 

Hie  parametric  decay  of  an  intense,  electromagnetic  wave  into  electro- 
static waves  and  scattered  electromagnetic  waves  in  fully  ionized  plasmas 
has  been  studied  extensively. 1-11  These  studies  have  wide  applications 
to  laser  fusion  experiments,  interactions  of  pulsar  radiation  with  a plasma, 
ionospheric  modification  experiments,  etc.  There  is  still  considerable  in- 
terest in  the  nonlinear  evolutionary  properties  of  these  waves  since  they 
ultimately  determine  the  importance  of  a particular  process  to  an  experi- 
ment. It  is  also  well  known  that  certain  nonlinearities  coepletely  change 
tiie  nature  of  a particular  mode;  e.g. , they  may  lead  to  explosive  instabil- 
ities. 

In  this  study,  we  examine  some  nonlinear  properties  associated  with 
the  panmetric  decay  of  an  intense,  coherent,  electromagnetic  wave  into 
an  electrostatic  wave,  its  second  harmonic  and  scattered  electromagnetic 
waves  in  a homogeneous  plasma.  For  the  case  of  scattering  off  ion  acoustic 
waves,  the  mode  coupling  processes  of  many  harmonics  must  generally  be  con- 
sidered. Nonetheless,  there  are  cases  in  which  a modulational  mode  and 
only  one  harmonic  can  exist  and  it  is  situations  of  this  kind  that  we  con- 
sider. 

We  find  that  general  conditions  exist  for  which  explosive  instabilities 
result,  and  that  this  phenomenon  may  occur  after  about  a linear  e-folding 
time. 


II.  GENERAL  FORMALISM 


Consider  the  propagation  of  a large -amplitude,  linearly  polarized, 
coherent  electromagnetic  pimp  wave, 

£ - 2 £ cos(£  »x  - Bt)  , 
o oo  o o 

in  a homogeneous  plasm.  In  the  absence  of  damping,  the  pair  (“0 »k0)  satis 
fies  the  usual  dispersion  relation. 


u a ■ k *c*  * o»  * , 
o o p * 

provided  that  the  electron  quiver  velocity  vQ  -(-efo/ntJo)  is  non-re  lativ- 
is tic. 

We  consider  the  case  in  which  the  ping)  wave  excites  a low-frequency 
aodulatianal  mode  (u,k)  and  its  harmonic  (2u,Zk)  , as  well  as  the  two 
lowest -order  high-frequency  sidebands  , (u±  ■ at  ± cu0,  t±  - £ ± kQ)  , 

(«2±  • 2w  ± <i>0,  * z£  ± £0)  , associated  with  each  wave.  The  effects  of 

pug)  depletion  are  neglected  and  it  is  assumed  that  all  waves  are  coherent. 

It  is  shown  below  that,  wider  these  conditions,  if  the  two  modes  satisfy 
certain  general  properties,  an  explosive  instability  results  and  the  waves 
become  wibounded  in  finite  time. 


First,  we  calculate  the  perturbed  current  densities,  since  they  are  the 
sources  that  determine  the  nonlinear  evolution  of  the  waves.  In  first  order, 
they  axe  as  follows: 


- at\  * £ PkXe(k,w) 


" ^SrT  tpl?e(k»“)  * ^Xi (*»“)]  . 
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where  o±  ■ ixjp2/4iru±  is  the  linear  plasma  conductivity  at  frequencies  oj±  , 
X(k,w)  - («p2^n0ka)  | dvk  • VvfQy{u>-it*v)  is  the  linear  susceptibility, 
and 

is  the  linear  field,  including  the  ponderomotive  contribution,  governing  the 
motion  of  the  electrons  at  frequency  u.  1 

In  second  order,  the  Vlasov  equations  for  the  perturbed  electron  density 
distributions  take  the  form: 

- ;[?k2)‘Vo  * VV-k  * Wa]  ■ » • 

-2i(«.-k-J)f®  - j[fg)-Vo  * VVk]  ■ » • 

- |[t0-Vk2)  * VV*0 

♦ vvf;a,‘)  * w*  * E-k-v2*  * y -v»  ] - ° • 

-iC“2»  - *21**2?  ■ jk-Va  ♦ W*a) 

L O 

^2,-i b 


(4) 


and 


*■?  ‘ffT 


tSa) 

(5b) 


(5c) 


(5d) 


(6a) 


(6b) 
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are  the  fields  that  determine  the  second-order  electron  dynamics.  18  The 
second -order,  perturbed  current  densities  are  due  to  the  electron  motion. 


They  can  be  reduced  to  the  following: 

if  - - * ^2) 

C7a) 

• - * ^2) 

(7b) 

*t<2)  ‘ 1 [S°  "k2>  * VkU*  * V “a”] 

(7c) 

^2)  - * ^ [W>  * eAa)] 

(7d) 

where 

"k1’  - *3“  V *'“>Pk  • 

(8a) 

»k2)  ’ WP  * Ml[s£  *.*•“>] ’Vk*] 

(8b) 

“a2)  ' I™-  [»e  'k>")pa2)  ' Hr  [sr  XgCktujJ  Pk>] 

(8c) 

are  the  perturbed  densities  corresponding  to  (w,k)  and  (2w,  Zk) 

in  lowest  order.* 

The  Fourier  - trans  formed  wave  equations  for  and  §+  take  the  form: 

- w*(l*X1Ck.M))Ejc  - ♦ 4iriwj^2^ 

(9) 

- • TT  Ve^’^k  * 47riwJ±2) 

(1C) 

where 

*f*  • f(ktc)J  - w±ae±]I  - k±5±c2  , 

(11) 

•See  Appendix. 


etc. 
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Due  to  linear  instability  and  nonlinear  coupling,  the  waves 
will  vary  on  a Hum  scale  roughly  of  the  order  of  a linear  e-folding  time.  The 
evolutionary  equations  for  the  wave  an|>litudes , obtained  in  the  usual  way , are 
as  follows:'*' 


aspt-iO^CZkJ-^OOM 


“R00 


X 


1 


ik^d+XiCk.uO)] 


(18a) 


I7 


«*p[-i(2cu^(k)  -u>^(2k))t] 


J(2oi.2k) 

(l+X-Czi.Zw)) 


t M^i*voi2ku0(l+xi(2k»2u))J 


where  * Real  (id)  and  the  fields  are  assumed  to  be  independent  of 
position.  17 


(18b) 


+Note  that  Eqs . (It)  art  be  modified  for  stimulated  Brillouin  scattering,  since 
in  that  case, 

A(«  ♦ i ,k)  - A(u»,k)  ♦ i ^ ^AGa.k)  ♦ ... 

-1  1 

does  not  converge  for  (3/3t)P^  « wP^  . This  is  because  M_  « 1/ (w-&w  ♦ ir_)  , 
and  o)  » dia  « kCs  for  ion  acoustic  modes,  so  that  M'1  cannot  be  expanded  in 
a Taylor  series  about  u . 
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Using  a linearized  version  of  Eq.(10),  J(w,k)  and  J(2w,2k)  can  be 
reduced  to  a term  proportional  to  P^P^*  ant*  Pk^  respectively,  where  the 
proportionality  constant  is  in  general  conplex.  Bqs.(18a),  (18b)  can  then 
be  transformed,  using  Eq.  (4) , into  a pair  of  equations  involving  and 
only,  'diich  take  the  form 

^St  “ Y^Pk  * ClP2kPk  » (19) 


(3t  ' VP2k  “ c2Pk 

t 

(20) 

respectively,  where 

i 

J 

(21a) 

' 

^Icl 

1 

r 

“R« 

Xe(lc,w)  ’ 

r « _ 

i 

J2 

(21b) 

2 

Xg  (2k,2tu)  » 

«R(2k) 

• 

J - JGu.kyp^* 

» 

(21c) 

J2  - J(2u,2k)/Pk2 

> 

(21d) 

the  small  frequency  mismatches  have  been  neglected,  and  terms  involving 
kv0*(M+1>  M_*)»kvQ  were  eliminated  by  using  the  linear  dispersion  relation, 

Eq.(17). 

A solution  to  Eqs . (19) 

,(20)  is  sufficient  to  determine  the  nonlinear 

evolution  of  all  wave  amplitudes,  provided  that  appropriate  initial  coreiiticms 
are  given.  11  They  apply  to  all  parametric  processes  for  which  both  A(cn,k)  * 

0 

and  A(2uf2k)  * 0. 
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III.  ANALYSIS 


Ne  now  examine  the  properties  of  Eqs . (19) , (20) . First,  we  transform 

i/ 

them  into  a dimensionless  polar  form.  Setting  T»yt,  A^»  [|cjC2|]  , 

and  Ajjj"  ciP2l/Y  * we  °*>ta^n 


( 3r  ' 1)Ak  * A2k\ 


(22) 


< h - «V 


(23) 


where  ♦ * argr(c^c2)  and  y ■ y(2k)/y(k).  In  the  polar  representation, 
Ak-  akel*k  # a^-  a^e^Zk  , . (22) , (23)  take  the  form 


JT  \ * ak(1  + a2k°°8  4,5 


(24) 


JFa2k  " *»2k  + •k2<»«(*-*)  . 


d V 

<t>  • — 3in(1>-<P)  - Za^sin  <J>  , 


(25) 

(26) 


where  ^ * Finally*  setting  x*a^e'2r  , y-a^e'^  , Eqs. (24) 

through  (26)  reduce  to  the  following: 


x - 2xy  ao8  <pe 


,rc 


(27) 


& - x <*>e(*-*)e<2^T  , 
y g|  ■ e(2‘^T[x  ainOp-*)-  2yV2(1'^T  exit*]  , 


(28) 

(29) 


which  uqply 

jp  xy  tin  <p  ■ x2e^2“^Tain  if> 

B^s.  (27) -(29)  are  now  in  a convenient  form  to  be  analyzed.  There  are  two 
general  cases  of  solutions  depending  upon  the  sign  of  ooe  <i>  . 

(i)  j at  » > 0 : 

We  divide  $ space  into  four  regions.  For  $ in  Quadrant  I , 
let  c ■ {$:(>  < ♦ < ♦)  and  c » < ♦ < ir}  , while  for  ♦ in  Quadrant  IV, 

set  c « {♦:*  < ♦ < 0}  and  c"«  (<>:0  < ♦ < <rnr}  . Let  c*-  c*  ir  and 

C * c + v . According  to  Eq.(29),  sgn(dWdr)  • ±1  for  + < <£- 

which  implies  that,  if  $ « c+  u c~  , then  $ changes  monotonically  until  it 


enters  c , where  it  is  trapped. 

For  ♦ e c , according  to  Eqs.(27),(28),  x and  y increase.  A lower 
bound  for  (x,y)  is  obtained  from  a solution  to 

jf  - 2xy  ooe  <|>  , 


^ - x ooe  i|i  , 


provided  that  f s 2.  This  corresponds  to  the  modulational  inodes  discussed  below. 
The  general  solution  to  Eqs. (31) , (32)  is  as  follows: 

^M^ton[X+  tan~l(y'/ M^2)]  for  M > 0 , 


jr’  ooth  X - |M|  ^etwh  X 1 % 

|M|  ^ooth  X - y'tinh  X J 


for  M < 0 , 


M - x - y2 


is  a constant  of  the  notion,  X » |M|  ^2  (T-x')ooe  $ , and  the  prise  refers  to 
to  an  initial  value  corresponding  to  the  tine  at  which  + enters  c . 


10 


According  to  Eqs. (33) ,(34) , x and  y become  unbounded  when 


[ir/2  - *m’V/Ml/2)yMl/2ooa  4 


T - T'  - T_  ■ 


— u — ■ — In 
2|M|  '2oo*4 


% 


y*-  mi' 

y'+  |M|  a 


for  M > 0 , 
for  M < 0 , 


(35) 


thus  providing  a conservative  estimate  of  the  time  at  which  x and  y become 
infinite.  For  $ < c*  , d4/dx  can  have  either  sign. 


(ii)  009  4 < 0: 

In  this  case,  let  c - (4:0  <♦<«/>>  • c“-  {4:4  < 4 < it), 
c ■ c ♦ n , and  c » c ♦ it  for  4 in  Quadrant  II,  and  let  c ■ {4:4  < 4 < 0)  , 
c‘-  (4:0  < 4 < 4+*}  , c*«c*ir  and  c+»c"+ir  for  * in  Quadrant  III. 

In  general,  for  $ e c u c+  , 4 changes  monotonically  as  before  until  it 

enters  c , where  it  is  "trapped",  and  in  general  oscillates. 

The  pair  (x,y)  , however,  may  become  explosively  unstable.  Assam  for 
definiteness  that  4 is  in  Quadrant  II.  Mien  * - w/2  < 4 < ir/2  , x and 
y both  increase  with  time.  Consider  the  region  4-4  <4  <4  for  some  4,  such  that 
4~ir/2  < 4 < ir/2  . A lower  bound  for  (x(T),y(x))  is  obtained  from  a solution 
to  the  equations 


jf  x - 2xy  <3oa  4 , (36) 

y " * <*>*  ♦ . (37) 

vdiich  is  described  by  Eqs. (33) -(35),  if  4 is  replaced  by  4.  The  pair 
(x,y)  therefore  becomes  explosively  unstable  if  4 remains  in  this  region 
for  a time  equal  to  t„  as  given  by  Eq. (35). 

If  4(o)  < c*  , d4/dr  can  have  either  sign. 
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Next,  we  calculate  J = J(u,k)  and  J^s  J(2w,2k) . Assuming  that 

^2k“  ^ * **  f3jld  thrt  *±  s±  » 'Aer®  tOTB  of  ar6ar  4 J^*«l 

are  neglected.  Also, 

* 4ire 

if  the  second  tern  in  Eq.(8b),  which  is  of  the  order 

a/iWiWiV 

8\ve/r/\“/  ’ 


"k  * s 4ro"  Xe  Pk2) 


is  neglected.  Finally,  comparing  the  two  terms  on  the  right-hand  side  of  Eq.  (16) , 
we  find  that  -die  second  tern  is  of  order 


/k\nv 

\ko)lc/ 


« 1 


of  the  first  and  can  therefore  also  be  discarded.  Since  these  results  apply 
to  , Eq.(16)  can  be  reduced  to  the  following: 


T . iek 


\ (.J*.  pk  e2k 

r , 

2 32nu^  \ “0/  *k 


(41a) 


(41b) 


where  ^ is  to  be  evaluated  at  (w,k)  . The  coupling  coefficients  take  the  form 


=1  * - T(w)^),x.,/CJIXil(1-<J/“i>,>  ♦ V'5' KT  * 2r/«)] 


(42a) 


°2  * -(»)(^)’x,‘/C2lxit<1-  * ix.t^  j^-  * r/u>]  . 


(42b) 


which  iaplies  that 


■ - ■-  ' - iiiigiiiifaMhji 


I 
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> » (w/kve)2  (43) 

for  an  explosive  instability,  where  r/u  was  neglected  coopered  to  u/k ve 
and,  since  Ix^l/x^  = «^a/«a  • Bq. (43)  is  a condition  which  can  be  readily 
satisfied. 

The  strength  of  the  nonlinearity  can  be  assessed  by  calculating  the 
initial  value  (x(o),y(o)).  According  to  the  definition 

x(o)  - |cjc2  P^/yI* 

, i /siW^oV is! pic 

5r)lv  r- 


where  the  inaginary  parts  of  and  c2  have  been  neglected  as  a rough 
appiuodnetion.  According  to  Eqs . (10) , (14a), 

'VE-I2  *(t=  *)“  • 

V 

Since  y/w^<1  , wa  * SI  ~ in  Gq.(39)  so  that 


14 


then,  for  an  initial  perturbation,  | E_/E  | a 2 x 10’*  , we  find  that  x(o)  * 1 , 


y(o)  a ✓|c1/c2l  a .3  , and  tp  a -19°  . Assuring  that  y'  - y(o)  in  Eq.(35), 
we  find  that  M • x(o)  - y*(o)  a .9  and  that  the  explosive  instability  occurs 
at  a time  $ Taa  1.4  linear  e-folding  times.  Since  (6/u)2«l  is  required 
for  these  modes,  we  note  that  harmonics  higher  than  the  second  do  not  satisfy 
the  linear  dispersion  relation  and  so  need  not  be  included  in  the  calculation. 


V.  DISCUSSION 

We  have  obtained  a general  formalism  describing  the  nonlinear  evolution 
of  the  system  (k,ui),  (k±  ,w±) , (2k,  2w),  Ck2±»<*>2±)  associated  with  the  para- 
metric decay  of  an  intense,  coherent  electromagnetic  wave.  This  system  was 
reduced  to  an  equivalent  two-wave  system  consisting  of  the  wave  fields  that 
determine  the  first-order  electron  dynamics,  and  . The  evolutionary 
equations  then  assumed  a standard  form  for  wave -wave  interactions  with  complex 
coupling  coefficients. 

The  polar  angle  , ■ arg (CjC2)  determines  the  essential  behavior  of 

the  system.  In  general,  the  system  oscillates  in  a "trapped  region"  of  t 
space.  If  t|>  lies  in  the  right-half  plane,  x and  y can  increase  simul- 
taneously and  an  explosive  instability  occurs. 

We  then  applied  our  results  to  a modulational  mode  and  found  that  it 
can  become  explosively  unstable  within  a time  of  the  order  of  a linear  e- 

I 

folding  time.  These  modes,  therefore,  can  grow  to  an  appreciable  size  of 
the  pump  (or  perhaps  of  any  other  modes  that  can  saturate  them)  within  this 
time. 
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Similarly, 


dv  P ,7  f 


- Va 


<Mc*v  f 


ie  f ”V0 

J (<ii4c»v*i6)  * 


ie  nok‘  1 3* 


pzkp-k  Ti^r'rssrxCk.uO 


(A.4) 


Bq*(A.l)  therefore  reduces  to  Eq. (8b).  Eq. (8c)  is  obtained  in  a similar  way. 
Next,  we  examine  . According  Xj  Eq. (5a), 

Jk2) ' m | £'Vo  * -a  Kf-k  * p-k£-VaJ 


e»k 

mTTiT 


dyjg.v-arid)^i8].  C2)^.7  f +p  {.7  £ *p  f , 
(uTTT^iS)  k Vo  P2k*V-k  P-k*VZk] 


--2f^2) 


(A.S) 


The  current  density  is  reduced  in  a similar  manner. 

The  side -band  contributions  ?±^  , 9X9  nore  arduous  to  obtain; 

hoMwer,  since  the  system  is  non-relativistic,  each  integral  is  easy  to  evalu- 
ate. It  is  necessary  to  have  expressions  for  the  linear  perturbed  electron 
density  distributions.  The  linear  Vlasov  equations  are 


- itv^fjW-  J£0*Vo  ’ 0 » CA.6) 

o 

- . a (8t.vvf0  . t0-7v^  ♦ VP*  > * 0 . c*.n 

o 

- i(»-5‘V)fk  - £ ?k*Vvf0  • 0 , (A.S) 

with  similar  expressions  for  f2k  * f2±  * 
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The  fourth  term 


e‘ 

S’ 


dv 


7 V V-t13 

(<*»±-Ic±*$*i<5) 


e^_ 

S’ 


dv 


-2± 


vflc  •£ 
lK±  c2± 


[(aj±-^±*v+ifi)  (w±4c±«v*i6)2J  *k 


:(1) 


!ii’p 

iau±  E2± 


(A. 12) 


in  the  lowest  order.  The  fifth  term 


e* 

S 


7 vyg1 

(«±4c±  •$♦!$) 


-g-f  J \ 

J |C»t-£±«v«l«)  (w±-Ic+*v+i5)2J 

i«  [ ^2±*Vq  * V7vf2k  * ^2k*V±k«  1 

" L K±*  J 


£j.  K^2t^2±*^ig)  * •$»!<) 

k.  (‘u±‘^+'^,ifi)2(h,2±"k2+*^ig^2 


<VB-k>  * 


(u±  -k±  • v+i6 ) 2 (U2±  -k2±  • v+ifi ) 


♦ ?(£t*2  p ^2±  ~^2jvfig)*k2tt(Ut~k±,^l‘i<s ) 

L (“± ~^± 'Y*W  (u^ ~^2±  'V+i& ) 2_ 

- tVo'Vzk  * Vtk ] 

o 

* * K7  ff-kC^VM-kllHVo-Vzk) 


(A.13) 
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